00100	COMMENT ⊗   VALID 00002 PAGES
00200	C REC  PAGE   DESCRIPTION
00300	C00001 00001
00400	C00002 00002	procedure J0J1(real mi reference real j0,j1)
00500	C00004 ENDMK
00600	C⊗;
     

00100	procedure J0J1(real mi; reference real j0,j1);
00200	begin
00300		real xs;
00400		fortran real procedure COS(real x);
00500		fortran real procedure SQRT(real x);
00600	
00700		if mi<3.0 then
00800		begin "L3"
00900		    xs←(mi/3)↑2;
01000		    j0←1-xs*(2.2499997+xs*(1.2656208+xs*(-0.3163866+
01100			xs*(0.0444479+xs*(-0.0039444+xs*0.00021)))));
01200		    j1←mi*(.5+xs*(-0.56249985+xs*(0.21093573+
01300			xs*(-0.03954289+xs*(0.00443319+xs*(-0.00031761+
01400			  xs*0.00001109))))));
01500		end "L3"
01600		else begin "GL3"
01700		    xs←3/mi;
01800		    j0←(1/sqrt(mi))*(0.79788456+xs*(-0.00000077+
01900			xs*(-0.0055274+xs*(-0.00009512+xs*(0.00137237+
02000			  xs*(-0.00072805+xs*0.00014476))))))*
02100			cos(mi-0.78539816+xs*(-0.04166397+
02200			  xs*(-0.00003954+xs*(0.00262573+
02300			    xs*(-0.00054125+xs*(-0.00029333+
02400			      xs*0.00013558))))));
02500		    j1←(1/sqrt(mi))*(0.79788456+xs*(0.00000156+
02600		      xs*(0.00017105+xs*(-0.00249511+xs*(0.00113653-
02700			xs*0.00020033)))))*
02800		      cos(mi-2.235619449+xs*(0.12499612+xs*(0.0000565+
02900			xs*(-0.00637879+xs*(0.00074348+xs*(0.00079824-
03000			  xs*0.00029166))))));
03100		end "GL3";
03200	end;